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We present a self consistent method based on cluster algorithms and Renormalization Group 
on the lattice to study critical systems numerically. We illustrate it by means of the 2D Ising 
model. We compute the critical exponents v and Tj and the renormalization group flow of the 
probability density function of the magnetization. The results, compared to the standard Monte 
Carlo Renormalization Group proposed by Swendsen |[TJ, are very accurate and the method works 
faster by a factor which grows monotonically with the lattice size. This allows to simulate larger 
lattices in reachable computational times. 
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1. The Method 

We consider a general spin model defined in a square lattice of lattice spacing a, and linear 
size L, with periodic boundary conditions which is defined by the Hamiltonian 

H = ^K a S a (1.1) 

a 

where each S a is a combination of the spin variables, for example Si = OjOj for first neighbors, 
S2=OiO m for second neighbors, and so on. 

In this article we illustrate the method by considering the 2D Ising model, including up to 
three even interactions (nearest-neighbor, second-neighbor, and four spin) and one odd interaction 
(a weak magnetic field). 

The RG theory allows the computation of the critical properties of a model. The critical 
exponents for example, can be obtained from the linearized RG transformation matrix T a p*, by 
computing its eigenvalues. T a p*, defined by 



T a t 



dK p (»- 1 ) 

can be obtained numerically from the coupled equations, 



(1.2) 



H* 



where, 



and 



a(s r w ) 



yj^E. i (13) 



S r (,l) Sp^-V } + ( S r W } ( Spl'-V } (1.4) 



^^>=-(S r (-)S a W> + (S r W)(S B W) (1.5) 

The critical exponents are obtained from the eigenvalues of T a p* in the standard way. For the 
Ising model they are given by the standard relations v = lns/lnAi 6 and rj =d + 2 — 21nAi°/ln&, 
where X\ e ^ is the largest even (odd) eigenvalue. 



2. The CARG method 

For the 2D Ising model the Hamiltonian reads 

%) = 7^I^i (2-D 

k » T <U> 

where / > describes a ferromagnetic system. In the thermodynamics limit this model has a 
second order phase transition at the Onsager's critical temperature T c = 2/ln(\ + \f(2)). 



2 



CARG Method 



Guillermo Palma 



2.1 Magnetic Susceptibility 



In order to compute the linearized RG transformation matrix T a p* from eqn. Ql.2[ ), the critical 
temperature should be used. One can obtain it by computing the susceptibility for different lattice 
sizes as it is shown in figure 1. Its maximum for each lattice size L defines the lattice shifted critical 
temperature T c (L). 




Figure 1: The maximum of the susceptibility % m [L\ is obtained using interpolation. For example, 
max(x m )[L = 108] = 2.2886, max(% m )[L = 64] = 2.3008, and max(x m )[L = 45] = 2.314. 



2.2 Finite Size Scaling 

Now we perform a finite size scaling analysis (FSS) of the numerical results obtained above to 
extract the infinite volume limit of the critical temperature. We use the standard analytical formula 
T C {L) = T c (oo) + a/L to fit the numerical values T C (L). The fit is shown in figure ^[ 




Figure 2: From the fit we obtain the numerical values T ( ■(°°) = 2.269, and a = 2.023. The convergence 
to T c (°°) is fairly fast, as can be seen from the following values for the lattice shifted critical temperature: 
T C (L = 108) = 2.2877, T C (L = 64) = 2.3006, T C (L = 45) = 2.3140. 
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The thermodynamic limit of the lattice shifted critical temperature T c (°°) agrees with the On- 
sager's critical value up to an error lest than 0.8 per one thousand. At this point we want to 
emphasize that a precise value for exact value of the critical temperature is no necessary to com- 
pute the eigenvalues of the linearized RG transformation matrix. In fact, it is only necessary to 



calculate the derivatives appealing in eqn ( |1.2| ) in a neighborhood of a "linear region", where they 
are essentially constant. 

2.3 Binder Cumulant 

One alternative method to compute the critical temperature of the infinite volume system is to 
use the Binder Cumulant 

/ x 1 ( M 4 ) 

^n = i-3^rf (2-2) 

which becomes independent of the lattice size at the Onsager's critical temperature. This behavior 
is shown in figure ||. 
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Figure 3: The intersection of these curves occurs at: T c (°°) — 2.26903 ±0.00059. This value differs in 1 
percent from Onsager's exact result. 

The result obtained with the Binder Cumulant is consistent with the above result obtained by 
a FSS analysis. 



3. Numerical Results 

3.1 Computing the critical exponents. 

To perform the simulation on the fundamental lattice we used the Wolff algorithm [g]. A total 
of 10 5 sweeps were used to thermalize the system and 10 6 configurations were used to compute 
thermal averages. Up to three even interactions (nearest-neighbor, second-neighbor and four spin 
or plaquette) and one odd interaction (magnetic field) were considered. We considered a lattice of 
lattice size L = 64 and performed two RG steps and used the renormalization parameter b = 2. The 
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results for the critical exponents v and 77 are displayed in the following table. Compared to the 
exact values, the numerical results agree remarkably well, with errors of one in one thousand. 



RG Step 


h e 


V 


V 


n 


1 


1.9586 


3.6856 


1.0311 


0.2362 


2 


1.9986 


3.6676 


1.0010 


0.2503 


Exact 


2 


3.6680 


1 


0.250 



Table 1: Critical exponents for L = 64 and T = 2.259. 

In order to compare our results with the ones obtained in reference [|3|], we considered a lattice 
of lattice size L = 108 and performed three RG steps using the renormalization parameter b = 3. 
The results are shown in the following table, where the label [S] stands for the result of ref. [||] and 
[PZ] for ours. 



RG Step 




V 


V 




1[5] 


2.852 


7.705 


1.048 


0.2828 


1[PZ] 


2.8635 


7.7062 


1.0443 


0.2825 


2[5] 


3.021 


7.828 


0.994 


0.2540 


2[pz] 


3.0027 


7.8269 


0.9992 


0.2542 


3[5] 


3.007 


7.831 


0.998 


0.2534 


3[pz] 


3.0013 


7.8361 


0.9996 


0.2521 


Exact 


3 


7.8452 


1 


0.250 



Table 2: The critical exponents obtained by the CARG method [PZ] and by Swendsen's method [S] for 
L = 108 are displayed. 

The numerical results for the critical exponents v and tj are very accurate and both results ([S] 
and [PZ]) agree remarkable well with their corresponding exact values. The main advantage of our 
method consists that it is faster than the method proposed by ^ by a factor which increases with 
the system size. 

4. RG-flowofthePDF 

4.1 Probability Density Function 

As it was pointed out in [0], the probability density function (PDF) of the magnetization corre- 
sponds to the partition function of an auxiliary theory, which equals the original theory plus a very 
small imaginary coupling. From RG theory one knows that partition functions are invariant under 
RG transformations. This is confirmed directly in the following figure. 

The first plot of figure |] shows the PDF for the fundamental lattice and for two coarse grained 
lattices. In the second plot the PDFs normalized to the first two moments are displayed. The three 
curves collapse onto one curve with high accuracy, showing the scale invariance of the PDF. 
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(m-<m>)/o 

Figure 4: L = 64 at T = T C {L) = 2.3008. A -» L = 64. A' -> L = 32. A" -> L = 16. 

5. Conclusions 

We have shown how to obtain very accurate values for the critical exponents and how to com- 
pute the RG flow of the PDF of the magnetization without the previous knowledge of the critical 
temperature. The method is self consistent and when compared to the original method proposed in 

our method is faster by a factor which grows linearly with the lattice size from 6 for L = 64 un- 
til a factor 10 for L = 162. The use of a cluster algorithm to simulate the fundamental Hamiltonian 
leads to this advantage, which allows to simulate larger lattices in reasonable computation times. 
A further application of the method in connection with universal fluctuations is reported in [S]. 
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